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ABSTRACT 

Planets have been observed in tight binary systems with separations less than 
20 AU. A likely formation scenario for such systems involves a dynamical capture, 
after which high relative inclinations are likely and may lead to Kozai oscillations. 
We numerically investigate the fate of an initially coplanar double-planet system 
in a class of binaries with separation ranging between 12 — 20 AU. Dynamical 
integrations of representative four-body systems are performed, each including a 
hot Jupiter and a second planet on a wider orbit. We find that, although such 
systems can remain stable at low relative inclinations (< 40°), high relative incli- 
nations are likely to lead to instabilities. This can be avoided if the planets are 
placed in a Kozai-stable zone within which mutual gravitational perturbations 
can suppress the Kozai mechanism. We investigate the possibility of inducing 
Kozai oscillations in the inner orbit by a weak coupling mechanism between the 
planets in which the coplanarity is broken due to a differential nodal preces- 
sion. Propagating perturbations from the stellar companion through a planetary 
system in this manner can have dramatic effects on the dynamical evolution of 
planetary systems, especially in tight binaries and can offer a reasonable expla- 
nation for eccentricity trends among planets observed in binary systems. We find 
that inducing such oscillations into the orbit of a hot Jupiter is more likely in 
tight binaries and an upper limit can be set on the binary separation above which 
these oscillations are not observed. 

Subject headings: binaries: close - - celestial mechanics - - methods: n-body 
simulations — planetary systems 



1. Introduction 



Among all extrasolar planets discovered to date , 20% are in multiple stellar systems 
( IDesidera fc Barbieril 120071 : lEggenberger et al.ll2007bl ). Binary stellar systems with circum- 
stellar and circumbinary debris disks have been observed with an overall occurrence rate 
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higher than that for debris disks around single stars ([Trilling et al.ll2006l ). Numerical sim- 
ulations suggest that protoplanetary disks embedded in binary star systems should b e able 
to fo rm both terrestrial and giant planets (jBarbieri et al.l 120021 ; I Quint ana et al.l |2002| ; iBoss 
2006h . 



Most binary and multiple st ellar systems found to be harboring planets are wide, with 
separations larger than 100 AU (lEggenberger et al.ll2004l ; iMugrauer et al.l I2004J ) . However, 
several systems with separations as low as ~ 20 AU have been shown to harbor giant planets. 
These are the b inary systems HD196885, 7 Cephei, and Gliese 86 ; and the higher order sys - 



tem HD41004 flCorreia et al.l 120081 : lHatzes et al.l 120031 : Els et al.ll200ll : IZucker et all 120041 ) . 
There are two important things to note about these syste ms. First, since radial- ve locity 



surveys have always been in biased against close binaries (lEggenberger et al.ll2007bl ). this 
sample probably underrepresents the frequency of planets in such binaries. In particular, the 
lower limit for the separation of binaries that can harbor planets may be smaller than ~ 20 
AU. Second, such systems constitute a unique data set which can allow us to test theoretical 
models for planet formation and evolution. This is because the presence of a stellar com- 
panion at such proximity can greatly influence these processes. A stellar companion as far 
as ~ 50 AU can weaken the chances for the formation of a giant p l anet by stirring, heating 
or tr u ncating protoplanet ary disks (jJang-Condelll 120071 : iBossI 120061 ; iMayer et al.ll2005l ; iKley 
200ll : iNelson et all hoOOal lbl) . Therefore, the fact that systems as tight as 20 AU are har- 
boring giant planets is puzzling and does seem to require further investigation. The most 
likely formation scenario for such systems may involve dynamical encounters in dense stellar 
systems where the binary companion is captured or moved inward from a wider orbit after 
the planet formation has taken place. A stellar companion can also affect the dynamical 
evolution of a planetary system, and therefore its long term survival, through secular per- 
turbations. In the case of a large relative inclination between the planetary orbit and that of 
the stellar companion, Kozai oscillations can take place where angular momentum exchange 
between the orbits leads t o large-amp li tude synchronous oscillations in the eccentricity and 



inclination of the planet (jKozail Il962l ; iHolman et al.l Il997l ; llnnanen et al.l Il997l ; iFord et al. 



20001 ; iTakeda &: Rasidl2005l ). This has the potential to disrupt the system if it results in 
close encounters, orbital crossings, or strong planet-planet interactions for systems with more 
than one planet. 

A renewed interest in planetary dynamics in general and in the stability of planets in 
binary systems particula rly has been triggered by the discovery of p l anets in binary and 



multiple stellar systems (|Dvorak et al. 



2003 



Barnes fc Raymond! 120041 ; iRaymond fc Barnes 



20051 : iRavmond. Barnes fc Kaiblhood iRivera fc Haghighipourl 120071 ). Looking into the dy- 



namical history of planetary systems is key to understanding the processes involved in their 
formation and subsequent evolution. In fact, this problem has received special interest from 
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dynamicists even before extrasolar planets were discovered, with many studies applying gen- 
eral results of the three body problem to the special case of a binary stellar system harb oring 



a planet (e.g, iDvorak fe Henrardlll993l ; llnnanen et al.l 119971 ; iHolman &: Wiegertlll999l . here- 



after HW99). By performing numerical simulations on planetary orbits in binary systems, 
HW99 derived empirical expressions for the maximum semi-major axis of the planet as a 
function of both binary mass fraction and eccentricity, above which the system would be- 
come unstable. Their study considered the stabilit y of binary system s harboring planets in 
the coplanar case only. Numerical simulations by iDavid et al.l (120031 ) of Earth-like planets 
in binary systems with different initial configurations led to the estimate that 50% of binary 
systems allow an Earth like planet to re main stable fo r 4.6 G yr. The stability of Earth-like 
planets in binaries was also studied by iFatuzzo et al.l (120061 ). With numerical simulations 
they calculated the survival time of an Earth-like planet orbiting a Sun-like star in the pres- 
ence of a stellar companion. They explored a range of binary config urations including re lative 
inclinations between the binary orbit and the orbit of the planet. lHaghighipourl (120061 ) per- 
formed numerical simulations to model the binary system 7 Cephei. This is a tight binary 
with separation 18.5 AU and an eccentricity of 0.36. A jovian planet wi th a minimum mas s 
of 1.7 M j was shown to be orbiting the primary of this system at 2.1 AU (IHatzes et al.ll2003l ). 
The question of the presence of an Earth-like planet within the habitable zone (HZ) of the 
primary was addressed and a range of semi-major axes and eccentricities of the binary as 
well as orbital inclination of the Jupiter-like planet was explored. This study concludes that, 
within timescales in the range 10 — 100 Myr, such a planet can exist on a stable orbit but 
not within the HZ of the primary. 

In this study, we investigate the stability of double-planet systems in tight binaries. 
The tightest stellar system in which a planet was ever claimed is HD188753. This is a 
hierarchical triple star system in which the primary is a Sun-like star with mass 1.06 M G . 
The binary companion is a double-star system with total m ass 1.63 M & orbiting the primary 
at separation ~ 12.3 AU and eccentricity 0.5 (jGriffinlll977l ). A hot Jupiter of mas s 1.14 Mj 
was claimed to be orbiting the primary star with a period of 3.35 d (lKonackill2005l ). Because 
of the proximity of the companion, the presence of a planet in this system is suggestive of a 
dynamical history in which the binary companion was either moved inward to a tighter orbit 
or captured after the claimed plan e t had already formed and reached its curr ent position 



(IPortegies Zwart fc McMillan! 120051 : IPfahll 12005 



Boss 



2006; Ijang-Condellll2007h . Although 



the observation of this planet was refuted by lEggenberger et al.l (l2007a[ ) , we find this system 
to be a good representative of the lower limit (in separation) of a class of tight binary systems 
harboring giant planets. We test in this study the capacity of such systems to harbor two 
planets including a hot Jupiter. This class of binaries with separation < 20 AU is of special 
interest since it is likely that such systems may have had active dynamical histories and 
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are expected to induce large perturbative forces to planetary systems with dramatic results. 
Studying such systems in the presence of more than one planet, can uncover the role played 
by binary companions in shaping planetary sys tems. This is espec ially important since multi- 
planet systems seem to be common in nature ( Fischer et al. 2003h. The dynamics of double - 
planet systems in wide binaries was studied analytically by iTakeda. Kita fc Rasiol (120081 ). 
where three different dynamical classes were identified as possible outcomes: (1) Decoupled 
systems in which planetary orbits experience independent Kozai cycles due to their weak 
mutual interaction compared to perturbations from the companion; (2) Weakly coupled 
systems in which mutual gravitational interactions become more dominant and relative nodal 
precession can cause mutual inclinations to grow and induce Kozai oscillations in the inner 
orbit; (3) Dynamically rigid systems in which the orbital elements of both planets oscillate in 
concert. These possible outcomes depend on the initial conditions and therefore the relative 
strength of the different perturbative forces. Here we expand on this study by looking into 
systems with closer companions and test the possibility of inducing eccentricity oscillations 
in a hot Jupiter orbit through a weakly coupled system. 

Since planets with orbital radii a < 0.1 AU are expected to be tidally circularized, the 
recent observations of a f ew hot Jupiters on considerably eccentric orbits (such as XO-3 , 
HAT-P-2, and HD185269 Jjohns-Krull et alil2008l : lLoeillet et allboosl : [johnson et al.ll2006h l 
does raise questions. Wh y are the orbits of these planets no t circularized when the orbits 
of most hot Jupiters are? iMatsumura. Takeda fc Rasiol (120081 ) show that close-in planets on 
eccentric orbits can be explained by constraining their tidal Q factors under the assumption 
that these orbits are in the process of being circularized. They also explore the possibility 
that these eccentricities are induced by a perturber on a wider orbit for the system GJ436 
and show that such a planet, if present, would cause a radial- velocity amplitude above the 
current detection limit, and therefore, would have been observed. 



Although such eccentricities can, in principle, be explained by the presence of a com- 
panion, whether stellar or sub-stellar, one has to keep in mind that a planet on a close-in 
orbit experiences strong GR forces. Therefore, a stellar companion as close as 12 — 20 AU 
cannot induce such eccentricities directly into a hot Jupiter orbit. We test here the possibil- 
ity of inducing such eccentricities to a hot Jupiter through a second planet in the presence 
of a binary companion. We also investigate the dynamics and stability criteria of a double- 
planet system orbiting the primary of a tight binary system with separation in the range 
ab = 12 — 20 AU. All systems we investigate in this study contain a hot Jupiter of mass 
lMj and separation a\ = 0.05 AU. The stability region for a second planet, placed on a 
wider orbit, that is initially circular and coplanar with a variable semi-major axis 02, was 
explored for two different binary configurations. In the first, the binary companion orbits at 
a>b = 12 AU with an eccentricity e& = 0.5, and in the second, a& = 20 AU and = 0. For 
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the tightest binary separation (a& = 12 AU), two different planetary masses rri2 = Mj and 
m 2 = M e were investigated. In this class of systems, planets experience mutual gravitational 
perturbations in addition to perturbations from the stellar companion. A hot Jupiter orbit 
is strongly dominated by GR forces. Because of the tightness of these systems, the interplay 
between these different forces is expected to have a dramatic effect on the planetary system. 

Throughout this paper, the subscripts 0, 1,2, and b will be used to refer to the primary 
star, the inner planet, the outer planet, and the binary companion, respectively. We use 
a semi-analytical approach to predict the stability region of the outer planet m 2 . Our 
numerical simulations incorporate a wide range of parameters including mutual inclinations 
of the orbits and perturbative forces to the Keplerian motion of each planet including planet- 
planet interactions and GR effects. Such numerical simulations can reveal to us the likelihood 
of a compact stellar system to maintain multiple planets and can be a useful guide for future 
observational efforts such as HARPS and the Kepler Mission. 

An overview of secular perturbations affecting this class of systems is presented in §2 . 
In §3 we discuss the analytical background for our predicted stability region for the outer 
planet in a stellar system resembling HD 188753. The numerical techniques used in this study 
are presented in §4, results and discussion in §5, and finally, §6 summarizes our conclusions 
and suggestions for future work. 



2. Secular Perturbations in Hierarchical Triple Systems 



Hierarchical triple systems can be studied using classical Hamiltonian perturbation tech- 
niques, which show that the secular evolution of the orbits takes th e form of eccentricity 
oscillations coupled with precession of the longitudes of pericenter f Brouwer &: Clemence 
196ll : lKozailll962l : lHarringtonlll968l : iRasio et al.lll997l : iFord et aljboooh . Such perturbations 
are highly dependent on the in itial relative inclination betw een the two orbits. The classical 
planetary perturbation theory (IBrouwer &: Clemencdll96ll ) applies to the case of low relative 
inclination orbits and low eccentricities for planets around a central star and is applicable to 
all orders in the ratio of the semi-major axes a, where a = ax/0,2, where a-\ and represen t 
the semi-major axes of the inner and the outer orbits respectively (e.g, lRasiolll994l . |l995). 



he s ecular evolution of a triple system with large relative inclination was derived by 



Kozai 



(119621 ) for the limit of small a using the quadrupole approximation of secular perturbation 



theory. In the following we discuss and calculate precession rates associated with these per- 
turbations in addition to GR forces as they apply to the type of systems investigated in this 
study (see §1). 
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2.1. The Low Inclination regime 

Mutual gravitational perturbations among planets can lead to a precession of the peri- 
astra and eccentricity oscillations of both planets. In the case of low eccentricities and low 
relative inclinations the periods of these oscillations, can be calculated using the formulation 
of the Laplace-Lagrange secular theory in which the disturbing potentials of the planets are 
simplified by eliminating short-period terms ( dependent on the mean long itudes). Following 



the classical planetary perturbation theory ( iBrouwer fc Clemencd ( 119611 ). see also Murray 



& Dermott 1999) for a two-planet system, the disturbing functions for the two planets, to 
second order in eccentricity, are 



n\a\ 



2^n e i + Ai 2 eie 2 cos(^i - u 2 ) 



(i; 



and 



n 2 a\ 



-A 22 e\ + A 2 ie 1 e 2 cos(co>i - uj 2 ) 



(2) 



where n, a, and u, refer to the mean motions, semi-major axes, and longitudes of pericenter 
of the planets, respectively. The matrix elements 



A 



ii 



1 m 2 

-ni ol b\L\a) 

4 m + mi 61 



(3) 



A 22 = \n 2 



mi 



4 m + m 2 61 



and 



An A 2 i 5 1 2 

• 4 ° = ^ = ^"-I a(1 -8 a) 



are functions of the planetary masses and semi- major axes, where a = ai/a 2 , and 



b { l(a) = - 

7T 



y 3/2 



2- 



cos^ dip 



o (1 -2a cos ip + a 2 ) 3 / 2 



(4) 
(5) 

(6) 



The eigenvalues of the matrix A represent the characteristic frequencies of the system 
with which the planetary orbital elements will evolve. These eigenfrequencies are given by 



9 



and 



1 
2 

1 
2 



(An + A 22 ) + J {An - A 22 f + AAqAhA 



22 



(7) 



(An - A 22 ) + J (A n - A 22 f + AAlAuA 



22 
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Pericenter precession rates are ,in general, non-li near combinations of these ei genvalues and 
can be estimated for different planetary systems (ITakeda. Kita Rasidl2008l ). The periods 
of oscillations can then be calculated from these rates. In this study, we use the period of 
orbital eccentricity oscillations given by P e - 



2- 



2.2. The High Inclination regime 



Kozail (119621 ) showed that in a hierarchical system where the ratio of semi-major axes, 
a, is sufficiently small and a large relative inclination (> 40°) between the orbits exists, 
exchange of angular momentum between the orbits results in simple periodic oscillations in 
both the inner eccentricity and the mutual inclination of the orbits i such that they are 
coupled by the integral of motion (1 — e 2 ) cos 2 2. Under these co nditions the amplitu de of 
eccentricity oscillations depends only on the relative inclination % (jlnnanen et al.l 119971 ) and 
is given at quadrupole order by, 



1 cos 2 % 

3 



(9) 



T he period of oscillati ons is a function of a, the masses, and the eccentricity of the outer 
orbit kiseleva et alJll998h 



(10) 



where P 2 is the orbital period of the planet. 



2.3. GR Forces 

For close-in planets, general relativistic effects become do minant and can cause the peri- 



astron of the planet's orbit to precess on very short timescales (jEggleton fc Kiseleva-Eggleton 



20011), 

5/2 / \ 3/2 

where a is the semi-major axis of the planet and mo is the mass of the central star. This 
precession can lead to the suppression of Kozai oscillations by averaging out the torques 
responsible for it if they act on shorter timescales. 
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2.4. Comparison of Timescales 

We compare in Figure 1 the precession rates due to all three perturbations acting on 
both planets in the case of a double-planet system in a binary with separation a& = 12 AU 
and eccentricity e& = 0.5. The planets, with masses m,\ = m 2 = Mj, are placed initially on 
circular and coplanar orbits with a\ = 0.05 AU. 

The dynamics of the outer planet are expected to be dominated by the Kozai effect as 
long as the semi-major axis of its orbit is larger than a 2 ~ 0.3 AU. Eccentricity oscillations 
with large amplitudes can cause this planet to experience close encounters either with the 
companion or with the hot Jupiter, leading to instabilities. Below a 2 « 0.3 AU, mutual 
gravitational interactions among the planets dominate, and Kozai oscillations are gradually 
suppressed. Therefore, a Kozai-stable zone exists below this point, within which this second 
planet can maintain stability even when relative initial inclinations are larger than the Kozai 
limit ii n it ~ 40°. The inner orbit at a\ = 0.05 AU is strongly affected by GR forces such that 
lugr > ^Koz and Kozai oscillations are suppressed in this planet. On the other hand, mutual 
planetary interactions can dominate the dynamics of the system as the outer planetary orbit 
gets tighter than a 2 < 0.2 AU and periastron precession rates associated with this force u pp 
exceed those associated with GR forces Cogr- This should be revealed in the form of mild 
eccentricity oscillations in both orbits. 

We look into the possibility of inducing Kozai oscillations in the inner orbit by prop- 
agating perturbations from the companion through the outer planet in the system. Kozai 
oscillations can be induced in the inner orbit if a relative nodal precession occurs, splitting 
the planes of the planetary orbits and producing relative inclinations larger than the Kozai 
critical angle. For these oscillations to occur two conditions must be satisfied: (1) The outer 
planet's orbit must experience nodal precession causing periodic growth in the mutual incli- 
nation angle between the planets. The outer planet can experience nodal precession even if it 
is not experiencing Kozai oscillations sinc e this precession is a natural consequence of secular 



interactions with the stellar companion (ITakeda. Kita fc Rasiol 120081 ) ; (2) Precession rates 



associated with the induced Kozai oscillation of the inner orbit must exceed GR precession 
rates. Whether both conditions can be satisfied simultaneously depends on the choice of 
initial conditions. In our systems the inner planetary orbit is fixed at 0.05 AU. The outer 
planet must be placed in a region where it is close enough for the induced Kozai effect to 
compete with GR forces. We show in Figure 2 the timescales for precession in the inner 
orbit caused by GR forces (eq.[ll]) and those caused by the mutual planetary Kozai effect 
calculated from equation (10) by replacing m& with m 2 as the perturber and replacing m 2 
with mi. Two different perturber masses are tested m 2 = Mj and m 2 = M®. It shows that 
while an Earth-mass planet cannot compete with the strength of GR forces acting on such 
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Fig. 1. — Secular precession rates, calculated from the periods such that Co = 2ir/P, for all 
perturbations acting on the inner planet (a), and on the outer planet (b) in a double-planet 
system with m\ = m 2 = Mj, and a\ = 0.05 AU. The binary parameters are a& = 12 AU, 
th = 0.5. Perturbations include planet-planet interactions (p-p), GR effects, and the Kozai 
secular perturbation. 
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Fig. 2. — The rate of precession in the inner orbit due to the GR effect compared to that 
associated with the mutual planetary Kozai effect for two different perturber masses m 2 = 
Mj and rri2 = M®. System parameters are the same as in Figure 1. The shown forces are: 
the GR effect (solid), the mutual planetary Kozai effect for m 2 = Mj (long-dashed), and for 
m 2 = M e (short-dashed). For equal mass planets, these two timescales become comparable 
at a 2 ~ 0.8 AU beyond which the induced Kozai cycles disappear. In the case of an Earth- 
mass planet in the outer orbit, mutual planetary Kozai perturbations become too weak to 
overcome GR forces. 
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a tight inner orbit, a Jupiter-mass planet can produce this effect when a 2 < 0.8 AU. 

3. Analytical Predictions for Stability of a Double-Planet System in a tight 

Binary 

We investigate the stability criteria of a double-planet system in a binary with separation 
a b < 20 AU. The inner planet is a hot Jupiter with semi-major axis a\ = 0.05 AU and mass 
mi = Mj. An approximate stability region for the outer planet is derived analytically 
within the two dimensional space [^4,02] where a 2 is its semi-major axis and ii n u is the 
initial relative inclination between the common planetary orbit and that of the binary. This 
stability region is derived analytically for two binary separations: a b = 12 AU and a b = 20 
AU. In the case a& = 12 AU, the stability region is compared for two different planetary 
masses m 2 = Mj and m 2 = M e in the outer orbit. For a given value of iinit, a planetary 
orbit is expected to be stable if a 2 lies between two limits, a 2i min and a 2jUlllx . The lower 
limit, a 2jiain , is defined as the minimum semi-major axis below which strong perturbations 
with the inner orbit can lead to instabilities, while the upper limit, c^max, is the maximum 
semi-major axis of the planet beyond which the planet enters a region of instability due to 
the presence of the stellar companion. 

3.1. Upper Limit 

The upper limit for the semi-major axis of a planet in a binary system was determined 
numerically by HW99 for initially circular planetary orbits which are coplanar with the 
binary orbit. This was done for a range of binary mass fractions fib = m b/( m o + m b) (where 
mo is the mass of the primary and m& is the mass of the perturber) and eccentricities e b of 
the binary. For each set of values [fib , e&], the binary system was numerically integrated 
with test particles representing the planet, and a critical semi-major axis a c was determined 
beyond which a planetary orbit may become unstable. An analytical expression for a c was 
derived from these numerical results such that, 

a c « (0.464 - 0.380/i 6 - 0.631e 6 + 0.586e 6/ u 6 + 0.150e^ - 0.198e^ 6 )a 6 (12) 

where a b is the semi-major axis of the binary. 

Applying equation (12) to the outer planet in the binary systems we consider gives a value of 
a c 1.23 AU for the stellar system HD188753, and a c R3 5.48 AU for the equal mass binary 
with separation 20 AU, both in the coplanar case. When large relative inclinations (i > 40°) 
are present, we estimate for an elliptical orbit the upper limit on the planet's semi-major 
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axis using, 

02,max ~ T-i . 7 (13) 

^ -L ~T &max) 

where e max is the maximum eccentricity gained during Kozai oscillations (eq. [1]), and as- 
suming the value of a c calculated from equation (12) is the distance from the primary at 
apastron. 



3.2. Lower Limit 



The stability of three-body systems was discussed analytically in the literature in more 
than one context. Some aimed at the stability limits for a system composed of a central 

Gladman] 1993, hereafter G93 ), while others were mainly 



star orbited by two planets (e.g 

focused on triple stellar systems (e.g. iMardling fc Aarsethll200ll . hereafter MA01). 
G93 provides stability criteria for a system of one star with two planets. He derives a 
minimum separation between the two planets "<5" below which the system cannot maintain 
Hill stability. A system is considered Hill stable if the planets cannot experience close 
encounters at any time. For systems in which the central star m is much more massive 
than both planets, m x and m 2 (/ii = mi/m << 1 and fa = m2/m << 1), and for initially 
circular and nearly coplanar orbits, the minimum separation 5, expressed in units such that 
the semi-major axis of the inner planet is unity, is given to lowest order in the masses as 



5 ~ 2.40(//i + (Mi) 



1/3 



(14) 



For equal mass planets this reduces to 5 ~ 3/i 1 / 3 . For non-circular orbits, but small eccen- 
tricities (e < /i 1 / 3 ) and equal-mass planets, this becomes, 



+ 9^/3. 



(15) 



This reduces to the previous equal-mass case for zero eccentricities. 

G93 also derived the Hill stability criterion for the case of equal mass-planets and equal and 
arbitrarily large eccentricities {e\ = e<i = e), giving 




1/3 > 



(16) 



This expression is valid only for eccentricities larger than e ~ /W 3 , in which case the mass 
terms become negligible compared to the large eccentricity terms. Therefore, at large eccen- 
tricities 5 becomes independent of planetary masses. 
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One can apply this stability criterion to the two planets in our systems. Eccentric orbits 
are expected if the binary orbit is highly inclined relative to the common planetary plane. In 
this case, 6 is calculated from equation (16) by substituting e max from equation (1) for the 
eccentricity. This sets a lower limit on the minimum separation allowed between the planets, 
since in some cases, the outer planetary orbit may experience Kozai oscillations while the 
inner planet is shielded from this effect due to stronger perturbations (§2.4). To calculate 5 
for systems with relative inclinations less than 40°, assuming e\ = e 2 ~ zero, we use equation 
(15) in the case of equal mass planets (m% = m 2 = Mj) and equation (14) in the case of 
a second planet of Earth-mass {m 2 = M®). Once 5 is estimated, the minimum semi-major 
axis allowed for the second planet can be calculated as a 2min = + 5), where a\ is the 
semi-major axis of the inner planet. The difference in the value of a2, m in obtained when 
applying the G93 criterion to the two different values of m 2 used in this study, is ~ 0.01 
AU. This is small enough that in practice here we ignore it when considering our predicted 
stability region. 

The stability of a hierarchical triple configuration was discussed in MA01, who describe 
hierarchical coordinates for an inner binary with masses mo and mi, and orbit described 
by ei n and aj„ (the inner eccentricity and semi- major axis). A more distant object m 2 
orbits around the center of mass of the inner binary with eccentricity e out and periastron 
separation R per . They derive a critical value for R per /a in below which the triple configuration 
can become unstable and may experience the escape of one object, 

gg = 2 8 ( (i+go^Kj+^ y /5 ; (17) 



(i - e out y/* 

where q out = m 2 /(m + mi) is the outer mass ratio. This eq uation gives an uppe r limit 



which represents coplanar systems. Numerical simulations by iMardling fc Aarsethl (119991 ) 
show that non-coplanar systems tend to be more stable, and obtain from their simulations 
a reduction factor / = 1 — 0.3i/ir (where i is the relative inclination between the inner and 
outer orbit) to account for the increased stability of inclined orbits. 

This stability criterion can also be applied to the two planets around the primary star in 
HD 188753 to place a lower limit on the distance of periastron of the second planet. One can 
again substitute for e out the value of e max , obtaining a lower limit on the semi-major axis of 
the second planet. This gives 

d2,mm = 2.8 X dl — — . (18) 

I, -L "max ) 



The MA01 criterion is weakly dependent on the mass of the outer planet m 2 through the 
ratio q out = m 2 / (m + m x ) in equation (18). Therefore, the difference in the estimated value 
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of a2,min obtained when using the two different values of m-i {Mj\ M®) is very small and 
therefore is neglected. 



3.3. The Predicted Stability Region 

In Figure 3 we show the analytically predicted stability region for a double-planet system 
with mi — ^2 — Mj in HD188753 where a hot Jupiter at 0.05 AU represe nts the inner planet. 



The upper limit (a2. m ax) is constrained by the lHolman &: Wiegertl (119991 ) criteria from which 



a value of ~ 1.23 AU was obtained in the coplanar case. For the lower limit on a 2 , both 
the G93 and MA01 stability criterion were applied to the central star in HD188753 with 
two planetary orbits. The maximum eccentricity of the outer planet was calculated from 
equation (1) for every value of the initial inclination ii n u- For the range ii n a < 40°, the value 
Gmax = was adopted. Therefore, this range represents the expected stability criterion for 
initially circular orbits. The minimum semi-major axis allowed for the second planet derived 
from G93 is a2,mm ~ 0.065 AU, while the MA01 sets a more stringent limit of a^mm ~ 0.14 
AU. Applying the reduction factor (/ = 1 — 0.3i/7r) for % = 70°, reduces this value to ~ 0.12 
AU. Figure 3 also shows the Kozai Suppression Limit (KSL). The shaded area below this 
line is the predicted Kozai-stable zone within which a second planet is expected to survive 
Kozai oscillations even in the presence of high relative inclinations (\ n iu a i — 40°). 



4. Numerical Methods 



The dynamics of a double-planet system in a tight binary with separation in the range 
ab = 12 — 20 AU is studied. We look into the effect a second planet may have on the orbit 
of a hot Jupiter in the presence of a close-in binary companion. We also investigate the 
stability criteria of this class of systems for two different binary configurations. The first, 
taken as a representative of the tightest binary in this class, resembles the stellar system 
HD188753 with separation 12 AU and eccentricity 0.5. The second is a 20 AU equal-mass 
binary (mo = m;, = M & ) with zero eccentricity. In all our integrated systems a hot Jupiter 
of mass nil — Mj is placed on an initially circular orbit at 0.05 AU. We consider two cases 
for a second planet: a Jupiter mass planet (equal mass planets), and an Earth mass planet. 
We examine the possibility of a habitable planet in HD188753 in the light of these results. 
For integrations resembling HD188753, the semi-major axis of the second planet was taken in 
the range 0.055 AU < ci2 < 1.8 AU in bins of 0.1 or 0.2 AU in all regions excluding the closest 
to the inner orbit where we examine numerically the cases a 2 = 0.055, 0.06, 0.063, 0.066 and 
0.1 AU. In the case of the 20 AU binary, a 2 was varied between the values 0.06 AU and 6.0 
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Fig. 3. — The analytically predicted stability region for a second planet in HD 188753 
with a hot Jupiter at 0.05 AU. The shaded region represents the most stringent limits 
expected for a second pl anet of mass m 2 = Mj. The upper limit is constrained by the 
Holman fc Wiegert Jl999h cr iteria (short dashed long dashed) . The low er limit is constrained 
by the iMardling fc Aarsethl (120011 ) (dashed), and the iGladmanl (119931 ) (dotted (large e) and 
dotted dashed (small e)) criteria. The Kozai Suppression limit is shown in long dashed. 
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AU in bins of 1 or 0.5 AU. In all of our simulations, the initial planetary orbits are circular 
and coplanar. The common planetary orbital plane is chosen initially to make an angle with 
the binary orbit i init , which was varied between 0° and 70° in steps of 10° in most cases. 

We perform a total of 5460 simulations representing 273 different initial configurations. 
For each configuration, many realizations were constructed with different initial relative 
phases between the orbits, which were chosen randomly. The number of realizations varied 
between 5 and 100 depending on the configuration in hand. 



Integratio ns were pe r 



orme d using a modified version of the integration package MER- 



CURY6.2 (I Chambers! (119991 )). The code was modified to include GR effe cts as an additional 



force. We adopted for this force the standard textbook \ term (e.g, iMisner et al.1 Il973 
eq. 25.16). Since our integrations are expected to involve close encounters, we use for all 
of our integrations the Bulirsch-Stoer integrator. We monitor conservation of total energy 
and angular momentum for the duration of the integration. Energy degeneration was always 
maintained to < 10~ 5 and fractional angular momentum change was always smaller than 
< 10~ 7 , excluding runs performed with the second planet's orbit close to the inner orbit 
(02 < 0.1 AU). In these cases, the fractional angular momentum change would fall in the 
range 10~ 6 — 10 -5 . 

Since direct integrations of such tight systems involve many close encounters, and are 
computationally expensive, and since we are not interested in following the dynamical evolu- 
tion of these systems after they become unstable, we optimize our computations by stopping 
them as soon as instability develops. We test the timescale for the onset of instability by 
integrating a system just outside the expected stability region. The system chosen for this 
test is one with a Jupiter mass planet placed on an orbit with semi-major axis 02 = 1.4 
AU and a relative initial inclination of 40°. We show the instability timescale histogram 
for this system in Figure 4, which includes 363 different realizations. We find that 99% of 
these simulations hit instability by 6.6 x 10 5 years. Based on these results we chose for our 
integrations the integration time 6.6 x 10 5 years to decide stability for a Jupiter mass planet. 
This timescale corresponds to more than 10 4 binary cycles. We also integrate systems with 
a Jupiter mass planet for as long as 6.6 Myr when this planet is near the Kozai suppression 
limit, since these systems tend to survive longer than their other counterparts. 



5. Results and Discussion 

We test the stability region within a tight binary for a double-planet system composed 
of a Jupiter mass planet at 0.05 AU and a second planet on a wider orbit. Two binary 
separations are investigated, 12 AU and 20 AU. The stability region is tested for a Jupiter 
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Fig. 4. — Instability timescale histogram for a Jupiter mass planet just outside the expected 
stability region, with a 2 = 1.4 AU and i init = 40° in a binary resembling HD188753. "N" 
is normalized to one and represents the cumulative fraction of systems that have already 
become unstable. We find that 99% of systems which end up unstable hit instability before 
6.6 x 10 5 yr. 
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mass planet in the outer orbit for both binary separations. In the case of 12 AU separation, 
the stability region for a potentially habitable planet is also tested where m 2 = M®. 

5.1. A 12 AU Binary: The Case of HD188753 

5.1.1. Equal-Mass Planets 

The stability region of the double-planet system described above in a binary system 
resembling HD188753 is investigated. Figure 5 summarizes the results of all our simulations 
of this system spanning the whole expected region of stability for a second planet of Jupiter 
mass. It shows that as long as the orbits are highly inclined (> 40°), Kozai oscillations 
will always lead to instabilities in this system. For small relative inclinations i init < 40°, 
the simulations give a lower limit for the semi-major axis of this planet in the range 0.063 
AU < a2,min < 0.066 AU, which is consistent with the G93 predicted value of 0.065 AU. 
The upper limit lies in the range 1.3 AU < a 2jmax < 1.4 AU for the coplanar case, larger 
than the value predicted by HW99. These results are obtained using full 4-body simulations 
which include planetary interactions with all other particles in the system and therefore, 
are expected to set more realistic limits on stability boundaries than do studies that treat 
planets as massless point particles. 

Our results also confirm the existence of a Kozai-stable zone within which stability is 
maintained even in the presence of large relative inclinations. This is seen in the region 
a 2 = 0.1 — 0.3 AU. At a 2 = 0.3 AU, we can see a transitional region between stability and 
instability as the system approaches the edge of the Kozai-stable zone. We compare in Figure 
6 the time evolution of both planets in two different cases. In the first the outer planet is 
within the Kozai-stable zone at a 2 = 0.2 AU. Although its orbit is highly inclined relative 
to the orbit of the companion, we see stable orbits with no sign of growth. The second case 
represents the outer planet outside the Kozai-stable zone at a 2 = 1.0 AU , where mutual 
gravitational perturbations become too weak to suppress Kozai oscillations. This case leads 
to the ejection of the planet within less than 80, 000 years. 

Figure 5 also shows a region of the explored space (0.35 < a 2 < 0.8 AU, and i ini > 40°) 
within which the planets are weakly coupled and therefore Kozai oscillations are induced in 
the inner orbit by interaction with the outer planet. Large eccentricities are gained during 
these oscillations which eventually bring the inner planet to the surface of the star. We 
discuss this effect in detail in 55.3. 
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Fig. 5. — The stability region for a second planet in HD 188753 of mass 1712 = mi = Mj. 
Empty triangles represent stable systems while filled triangles represent systems which ex- 
perience instabilities. All systems with the binary orbital plane highly inclined relative to 
the planetary plane are unstable, except for systems within the Kozai-stable zone. 
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Fig. 6. — Orbital elements of both planets (inner (dotted) and outer (solid)) as they evolve 
with time for the following systems: - Left panel: (m 1 = m 2 = Mj), ai = 0.05 AU, a 2 = 0.2 
AU (within the Kozai-stable zone), and ii n u = 70°. - Right panel: Same masses, a 2 = 1.0 
AU , and i init = 50°. The figure shows the distance at pericenter (r p ) and the eccentricity 
(e) of the planets. 



5.1.2. A Terrestrial Planet 



There has been great interest late l y in the question of habitability of exop l anets (e.g., 
David et alil2003l : Ijianghui et al.l [2OO5I : iFatuzzo et al.lhood : lHaghighipourl l20od . I2OO7L Al- 
though observations of Earth mass planets within the Habitable Zones (HZ) of Sun-like stars 
may not become reality in the very near future, numerical simulations to understand the 
dynamics and stability limits of these planets in different environments is necessary to guide 
future efforts. 

It has been shown that Earth-like planets can survive the migration of a gi ant planet 
( iMandell. Raymond fc Sigurdssonl 120071 : iRaymond. Mandell fc Sigurdsson 1 120061 ) . There- 
fore, we look into the stability of a double-planet system with a habitable planet and a 
hot Jupiter in HD188753. The central star in this system with mass 1.06 M ffi is of G-type 
with a habitable zone (HZ) between 0.95 AU and 1.37 AU feasting. Whitmire. fc Reynolds 
19931 ). Therefore, our predicted stability region (Fig. 5) does leave space for a habitable 
planet in this system, but only for small values of the initial relative inclinations of the orbits. 

Figure 7 shows the results of our simulations performed for a second planet of mass 
m,2 = M© in HD188753. Since an Earth mass planet is three orders of magnitude less 
massive than a Jupiter mass, perturbations from the companion are expected to become 
less dominant leaving more space for stability. Therefore, a maximum integration time of 
6.6 x 10 6 yr was used for all systems shown to account for any delayed instabilities. Our results 
show that the upper limit for stability of an Earth mass planet at low relative inclinations 
{Unit < 40°) is larger than that for a Jupiter mass planet and lies in the range 1.4 AU 
< (o-2)maiE < 1-5 AU, confirming the dependence of perturbations, and therefore stability 
limits, on the planetary-masses. The HZ of the central star in HD 188753 lies completely 
within the stable region at low inclinations. On the other hand, the Kozai-stable zone is 
located outside the HZ, leaving no space for a habitable planet when high relative inclinations 
are present. At the lower limit of a 2 , although an Earth mass planet can survive down to a 
region between 0.06 AU and 0.063 AU such a planet will be extremely hot and inhabitable. 

In addition to dynamical stability, an important condition for habitability is climate 
stability which can be comprimised if the planet's orbit becomes eccentric. This has the 
potential to drive the planet's orbit outside the HZ at apastron or/and periastron, or may 
lead to extreme heat variations on the surface of the planet. We show in Figure 8 the orbital 
elements of an Earth-mass planet in this system with low relative inclination and initially 
on a circular orbit at a 2 ~ 1.0 AU. It experiences mild eccentricity oscillations with an 
amplitude of e ~ 0.1, with no sign of growth. Still, habitability of this planet can become 
questionable since such an eccentricity can drive the planet near the edges of the HZ at 
periastron and apastron. Therefore, the presence of a perturber to the Keplerian motion of 
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Fig. 7. — The results of all numerical simulations performed for a system similar to HD 188753 
with a second planet of mass = 1M© . Empty triangles represent stable systems while filled 
triangles represent unstable systems. The stable region includes the HZ of this star, but only 
for low relative inclinations among the orbits. 
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Fig. 8. — Orbital elements of the outer planet in a double planet system in HD188753 with 
the outer planet in the Habitable Zone of the primary. System elements are : a\ = 0.05 AU; 
a 2 = 1.0 AU; m x = Mj\ m 2 = M e and i init = 0° . The figure shows from top to bottom: 
semi-major axis, eccentricity and argument of pericenter of the outer planet. 
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a potentially habitable planet can alter the positions of the outer and the inner edges of the 
HZ through their effect on the orbit. In this case, the outer edge of the HZ is expected to 
move inward to a 2 ~ 1.24 AU from the primary while the inner edge move from a 2 ~ 0.95 
AU to ~ 1.0 AU, altering the width of the HZ by ~ 40%. 



5.2. A 20 AU Equal mass Binary 



It has been shown by I Jang-Condelll (120071 ) that a stellar system as tight as HD188753 is 



unlikely to support a disk massive enough to form a Jovian planet in situ. Therefore, the pres- 
ence of a planet in such a system would be indicative of a dynamical history. On the contrary, 
the truncated disk around a binary with slightly larg er separation such as 7 Cephei (sep- 



arati on 20 AU and eccentricity 0.4) has been shown (IJang-Condell. Mugrauer. fc Schmidt 



20081 ) to have sufficient mass to form such a planet. Such wider binaries may host giant 
planets and still have no active dynamical histories, and therefore, maintain dynamical sta- 
bility with multiple planets which are coplanar with the binary orbit. Therefore, a binary 
with separation 20 AU, represents an upper limit in this study to a class of binaries with 
low separation in which the presence of a hot Jupiter is questionable. 

We show in this section, our results for the stability region of a double planet system in 
a binary with two stars of solar mass orbiting at 20 AU. The planetary system used in this 
test is the same system used in §5.1.1 in which two coplanar, equal mass planets of Jupiter 
size are place on circular orbits. The inner planet is a hot Jupiter at 0.05 AU, while the semi- 
major axis of the second planet, a 2 , is varied along the expected range. The stability region 
is also tested here in the two dimensional space [ii n u, 0,2] where Zj n # is the initial relative 
inclination between the binary orbit and the common planetary plane. We apply the same 
criteria as in §3 to predict the stability region for the outer planet when the inner orbit is 
fixed. These results are shown in Figure 9 where one can see that the stability region for this 
system is similar to that of HD 188753 in its structure but varies in scale. The larger binary 
separation has allowed more space for the outer planet to survive. In some cases the outer 
planet can experience Kozai oscillations (with moderate eccentricities) for as long as 1 Myr 
without disrupting the system as long as it does not reach the upper limit for stability at 
apastron. The Kozai stable zone expands with the binary separation in general since mutual 
planetary interactions become more dominant relative to the Kozai effect. In this system 
the edge of the Kozai stable zone is at a 2 = 0.5 AU. We also see that the two planets are 
placed in the "weakly coupled" dynamical class in a small region of the shown parameter 
space (0.5 > a 2 < 0.8 AU and %i n it = 50°) within which the inner planet experiences Kozai 
oscillations induced by the outer planet. The width of this region is narrower than that 
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observed in the case a& = 12 AU (see §5.3 for a more detailed discussion). 

5.3. Induced Eccentricity Oscillations in a hot Jupiter 

We test the possibility of inducing Kozai oscillations in the hot Jupiter orbit in our 
systems by propagating perturbations from the stellar companion through a second planet 
of equal mass. We find that inducing large eccentricities in the hot Jupiter orbit is more likely 
in binaries with smaller separation. In the case of a 12 AU binary separation this effect is 
observed within the region of the explored phase space defined by (0.35 < a 2 > 0.8 AU, and 
Uni > 40°). This behavior is not observed within the Kozai-stable zone (a 2 < 0.3) AU, and 
shows its strongest behavior in the region (0.4 AU < a 2 < 0.5 AU). As the distance between 
the planets increases, this effect gets weaker until it disappears beyond a 2 > 0.8 AU. Figure 
10 shows the time evolution of the orbital elements of a system representative of this behavior 
with the outer planet at a 2 = 0.5 AU and i init = 50°. The outer planet is experiencing Kozai 
oscillations induced by the stellar companion. The nodal precession of this planet creates 
a split in the planetary orbits and the mutual inclination between the planets, I\ 2 , grows 
and oscillates with the same period as fl 2 . The amplitude of these oscillations is twice the 
common initial inclination relative to the binary plane iinit- Angular momentum exchanged 
between the planetary orbits, propagating the angular momentum transfer from the binary 
leads to mutual planetary Kozai oscillations. Since the angular momentum of the perturber 
(m 2 in this case) is changing, the eccentricity oscillations induced in the inner orbit lead to 
chaotic growth. As the planet approaches the central star, it should experience strong tidal 
dissipation and may be circularized into a tighter orbit. Our simulations do not include tidal 
effects and therefore, we cannot predict the final outcome of these events. In general, tidal 
torques are not significant in the evolution of a circular orbit at ~ 0.05 AU where GR is the 
most dominant force. Once a planet gains significant eccentricities driving it very close to 
the surface of the star at periastron, tidal forces may become a dominant factor but do act 
on longer timescales when compared to Kozai precession timescales. 

To demonstrate the interplay between the different forces acting on the inner planet, we 
show in Figure 11 the relative inclination between the planetary orbits and the eccentricities 
of both planets in a system similar to that shown in Figure 10, but once with the GR forces 
removed from the calculation (top two panels) and second using an Earth-mass planet in 
the outer orbit (lower two panels). The former shows the faster growth of eccentricity when 
compared to Figure 10 and therefore demonstrates the goal played by GR forces in stabilizing 
the inner orbit against other perturbations, while the later confirms the dependency of this 
mutual planetary Kozai effect on the planetary masses. 
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Fig. 9. — The stability region for a double planet system in a binary with parameters: Two 
stars of Solar mass, a& = 20 AU, and = 0. Planetary masses are m 2 = m\ = Mj initially 
on circular coplanar orbits. The inner planet is a hot Jupiter at a 2 = 0.05 AU. Empty 
triangles represent stable systems, while filled triangles represent the unstable systems. 
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Fig. 10. — A case representative of the systems which experience induced eccentricity in the 
inner orbit, showing orbital elements of both planets, inner (dotted line) and outer (solid 
line). The initial system parameters are: mi = 1712 = Mj, a\ = 0.05 AU, a<i = 0.5 AU, 
t\ = e 2 = 0, and i init = 50°. The outer planet is experiencing Kozai oscillations induced 
by the stellar companion. Due to fast nodal precession of the outer planet, the mutual 
inclination between the planets Iu grows, and oscillates with the same period as 1^2- The 
amplitude of these oscillations grows to a value twice the common initial inclination relative 
to the binary plane, i init , causing the inner planet's eccentricity to grow chaotically. 
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Fig. 11. — Eccentricity of both planets (outer in (solid) and inner in (dotted)) and mutual 
inclination between planets I±2 shown for the same system in Figure 10, but with GR forces 
absent in the top two panels and and m 2 = M e in the lower two panels. 
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Inducing such oscillations into the orbit of a hot Jupiter requires specific initial con- 
ditions to be satisfied (see §2.4). We find that the region in parameter space that allows 
these conditions to be satisfied is narrower for wider binaries. We show in Figure 12 the 
range in semi- major axis of the outer planet that will allow induced Kozai oscillations in 
the inner orbit as a function of binary separation. This is shown for an inner planet on an 
initially circular orbit at a\ = 0.05 AU, and a binary orbit with inclination ii n u = 50° relative 
to the common planetary plane. All binaries in this figure share the following parameters 
(m = m b = M Q , mi = m 2 = Mj , and e& = 0). These results show that tighter binaries are 
more likely to produce this effect and that a cut-off binary separation a& iCUt exists beyond 
which it is not observed. For a hot Jupiter of mass Mj and semi-major axis a 2 = 0.05 AU, 
this cut-off binary separation lies in the range 30 < a& iCUt < 40 AU. As the inner planet's 
orbit is made wider, mutual planetary interactions grow larger which should allow a wider 
binary to induce eccentricities through weak coupling between the planets into the inner 
orbit. 



6. Summary and Direction for Future Work 

The presence of a giant planet in a compact stellar system with separation < 20 AU is 
indicative of an active dynamical history in which the binary companion is either captured 
or moved inward to form a harder binary during an encounter after the planet has already 
formed. Such an event is likely to result in high relative inclinations between the orbits of 
the planet and the companion and may lead to Kozai oscillations in the planetary orbit /s. 
Since multi-planet systems seem to be common, it is likely that the primary in such a 
system may be harboring more than one planet at the time of its encounter with a stellar 
companion. In such a scenario, Kozai oscillations can easily lead to instabilities including 
collisional events and ejections. Therefore, such encounters which may occur in dense stellar 
systems when planetary systems are in their infancy, can be very important events in shaping 
the final outcome of planetary systems and leading to the large diversity observed in their 
characteristics. 

We investigate the dynamical evolution of a double-planet system in a class of tight 
binaries with separation in the range 12 — 20 AU assuming relative inclinations between the 
common planetary plane and that of the binary in the range — 70°. Our simulations reveal 
the large parameter space involved with a double-planet system in a binary such that the 
final outcome of a system is highly sensitive to the initial conditions. The stability criteria 
are tested for two systems with binary separations 12 and 20 AU. We find that if the orbit 
of the stellar companion is highly inclined (> 40°) relative to the common planetary plane, 
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Fig. 12. — The region in parameter space which allows Kozai oscillations to be induced in a 
hot Jupiter orbiting at 0.05 AU. On the x-axis, a& represents the variable binary separation. 
Binary parameters for all binaries are: itlq = rrib = M Q ,mi = vri2 = Mj,eb = 0, and 
Unit — 50°. The variable semi-major axis of the outer planet 02 is plotted on the y-axis. 
In order to overcome GR precession and induce Kozai oscillations in the inner orbit, the 
outer planet must lie below the line a 2 = 0.8 AU. The dashed line represents the limit of the 
Kozai-stable zone as a function of binary separation. Inside the dashed region between the 
two limits, both planets experience Kozai oscillations (outer planet induced by the binary 
and the inner planet induced by the outer planet). All circles represent results of numerical 
simulations. The large filled circles are systems in which the outer planet experiences Kozai 
oscillations and always lead to eccentric hot Jupiters. Small filled circles are systems in which 
the outer planet is shielded from the Kozai effect, but in less than 50% of the simulations 
we observe induced eccentricities to the hot Jupiter. Empty circles represent systems which 
did not show any sign of this effect. 
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the system is likely to experience instabilities unless the outer planet is situated in a narrow 
Kozai-stable zone within which gravitational perturbations with the inner planet are strong 
enough to suppress the Kozai mechanism. At inclinations smaller than the critical Kozai 
angle, we find that a binary as tight as 12 AU can remain stable and an Earth-like planet 
can survive within the HZ of the central star up to 6 million years. 

The eccentric hot Jupiter phenomenon is a puzzling observation since hot Jupiters are 
expected to be the result of a migration process during which circularization of the orbit 
occurs. We suggest that eccentricities can be pumped up into a hot Jupiter by inducing 
Kozai oscillations into its orbit through a weak coupling mechanism with a second planet 
in the presence of a close-in binary companion. A fine tuning of the initial conditions can 
produce this effect allowing perturbations from the stellar companion to propagate through 
the planetary system into the inner orbit. We find that tighter binaries are more likely to 
produce this effect, and that for a given hot Jupiter orbit, a cut-off binary separation can 
be found beyond which inducing Kozai oscillation into its orbit becomes impossible. The 
final outcome of such a scenario as the inner planet gains large eccentricities cannot be 
predicted accurately without the inclusion of tidal torques. These forces become important 
as the planet approaches the primary at periastron during cycles of large eccentricities. 
Circularization into a tighter orbit is one possibility that seems unlikely if the outer planet 
remains in the system undergoing nodal precession and inducing oscillations in the inner 
orbit. This is since the timescales associated with these oscillations are orders of magnitude 
shorter than tidal circularization timescales. If the outer planet reaches the edge of the Hill's 
sphere of the primary at apastron and is ejected from the system, the inner orbit may have 
enough time to experience tidal circularization during the life-time of the star. Explaining 
eccentric hot Jupiters by such a mechanism requires the inclusion of tidal circularization 
effects in more comprehensive studies of binaries harboring double-planet systems. Such 
studies may require numerical testing over longer timescales. This mechanism can also be 
used to explain the higher eccentricities among planets hosted by binary systems in general 
and the relativ e paucity of planets in bi naries with near-circular orbits (for periods greater 



than 40 days) flEggenberger et al.ll2004j ). This is because the wider the orbit of the inner 
planet the easier it is to induce eccentricities into it by weak coupling with another planet 
on a wider orbit. 

Tight binary systems with planets constitute a unique class of binaries which can allow 
us to test theoretical models of planetary formation and evolution. Therefore, more com- 
prehensive studies including multiple planets and stellar collisions are necessary for a better 
understanding of the forces that were responsible for shaping planetary systems with a large 
diversity in characteristics. 

Finally, we can benefit greatly from better observational constraints on these systems, espe- 
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cially orbital inclinations. 
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